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A promising method for measuring the total mass of the dark matter near a supermassive black 
hole at the center of the Galaxy based on observations of nonrelativistic precession of the orbits 
of fast SO stars together with constraints on the annihilation signal from the dark matter particles 
has been discussed. An analytical expression for the precession angle has been obtained under the 
assumption of a power-law profile of the dark matter density. In the near future, modern telescopes 
will be able to measure the precession of the orbits of SO stars or to obtain a strong bound on it. The 
mass of the dark matter necessary for the explanation of the observed excess of gamma radiation 
owing to the annihilation of the dark matter particles has been calculated with allowance for the 
Sommerfeld effect. 


Significant advances have been achieved in recent years 
in the observations of stars gravitationally connected to 
the supermassive black hole SgrA* at the center of the 
Galaxy. Several so-called SO stars, which move at very 
high velocities (> 10^ km/s) in almost elliptic orbits 
around a very compact smermassive object, are observed 
in the infrared range [ll-lj ■ Models alternative to the su¬ 
permassive black hole at the center of the Galaxy, e.g., 
a cluster of compact stars such as white dwarfs, neu¬ 
tron stars, or black holes with star masses, are almost 
excluded [3, 0 • The SgrA* at the center of the Galaxy 
is most probably a supermassive black hole, although for 
the final proof it is necessary to confirm the existence of 
the event horizon of this object. 

According to the measured parameters of the Kepler 
orbits of SO stars, the mass of supermassive black hole 
SgrA* is Mbh = (4.1 ±0.4) x 10^Mq [il-[3. Independent 
and currently most accurate values of the mass Mbh and 
spin (Kerr parameter) a of the SgrA* black hole are de¬ 
termined from the observations of quasiperiodic oscilla¬ 
tions with average periods of 11.5 and 19 min [3,0. Tl^ 
are Mbh = (4.2 ± 0.2) x IO^Mq and a = 0.65 ± 0.05 Q. 

At the center of the Galaxy, in addition to the super¬ 
massive black hole SgrA*, there are additional invisible 
sources of mass such as compact gas clouds, dim stars 
and their remnants, and a distributed mass in the form 
of the dark matter density peak. Constraints on the dark 
matter density at the center of the Galaxy based on pul¬ 
sar effects were discussed in naini- All this additional 
mass would result in the deviation of the total Newto¬ 
nian gravitational potential from the potential of a point 
mass of the black hole U = —GMBn/r. As a result, the 
orbits of SO stars gravitationally connected to the black 
hole would be unclosed and precess (see, e.g., [ 10 ). The 
openness of the orbit of the most studied SO-2 star will 
be measured in the next one or two years. Thus, the 
total mass of the dark matter within the orbit of this 
star with a characteristic radius of 0.005 pc will be de¬ 
termined. The nonrelativistic precession of orbits of fast 
SO stars under consideration, depending on the mass of 
the dark matter near the center of the Galaxy, can sig¬ 


nificantly exceed the corresponding relativistic precession 
(an effect such as the shift of the perihelion of Mercury 
and frame dragging). 

The existence of fast SO stars provides a unique possi¬ 
bility of reconstructing the gravitational potential and 
measuring the mass distribution at the center of the 
Galaxy by fitting their orbits. The authors of EH per¬ 
formed a detailed multiparametric fitting of the orbits of 
several SO stars and calculated the additional distributed 
mass with various exponents of the density profile. It 
was shown that the distributed mass within the orbit of 
the SO-2 star is no more than 3 — 4 % of the mass of the 
supermassive black hole. It is noteworthy that the ex¬ 
pected measurement of the nonrelativistic precession of 
the orbit of the SO-2 star will allow either improving the 
indicated bound on the distributed dark mass by two or 
three orders of magnitude or determining this dark mass. 
We discuss and develop a method for studying the dis¬ 
tribution of the dark matter at the center of the Galaxy 
by measuring the precession angle of orbits of SO stars. 
For a number of particular cases, numerical calculations 
of the precession angle of orbits of SO stars because of 
the extended mass distribution were performed [l^ - [l0 . 
We obtained general analytical formulas for the preces¬ 
sion of orbits of stars with a powerlaw profile of the dark 
matter; these formulas make it possible to easily deter¬ 
mine the additional distributed mass from the measured 
precession angle. 

An additional independent method for determining the 
distribution of the dark matter is the search for a possible 
annihilation signal from the center of the Galaxy. The 
explanation of the excess of a gamma signal with an en¬ 
ergy of ^ 1 TeV from the center of the Galaxy observed 
by the HESS telescope by gamma annihilation of dark 
matter particles with allowance for constraints on the 
dynamics of stars for the case of the power-law density 
profile of the dark matter with a spike and an exponent 
as a free parameter was analysed in 2]|. The possibility 
of constraints on annihilation based on the dynamics of 
stars or precession was also mentioned in [ 10 . We cal¬ 
culated (see Figs. [T] and [U the mass of the dark matter 
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FIG. 1. Shift angle of the apsis of the orbit of the star in 
one turn Scj) calculated by Eq. © versus the exponent of the 
power-law spectrum of the dark matter j3 in Eq. © for real¬ 
istic values of the mass fraction of the dark matter ^ within 
the orbit of the SO-2 star. The indicated region is excluded 
by the constraints caused by the annihilation of dark matter 
particles if the dark matter makes the main contribution to 
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(see [13, Sect. 15, Problem 3) 



Here, integration is performed with the trajectory of 
the particle in the form of an unperturbed elliptic or¬ 
bit r((/)) = p{l -k ecos0)“^, where e is the eccentricity 
of the ellipse, p = L?/{GM^Yi'm) = a(l — e^) is the pa¬ 
rameter of the orbit, a is the major semiaxis, and L is 
the conserved angular momentum of the star with the 
mass TO. The observed parameters of the Kepler orbit of 
the SO-2 star: the eccentricity e = 0.898 ± 0.0034, the 
radius of the pericenter Vp = a(l — e) = 0.585 mpc, and 
the radius of the apocenter = a(l + e) = 9.42 mpc. 
We note that, in the case of relativistic precession, the 
orbit would rotate in the direction of the rotation of the 
star, but Newtonian precession m occurs in the opposite 
direction, i.e., S4> < 0. 

We consider the power-law density profile of matter 
responsible for the correction SU to the potential of the 
black hole: 
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FIG. 2. Mass fraction of the dark matter ^ versus the expo¬ 
nent P in the density profile given by Eq. © at the precession 
angle Sp = 0.01. The indicated region is excluded by the con¬ 
straints caused by the annihilation of dark matter particles. 


necessary for the explanation of the excess of gamma ra¬ 
diation from the center of the Galaxy detected recently 
by the Fermi-LAT space gamma telescope [2^ [2^. In 
particular, we determined the dependence of the addi¬ 
tional mass both on the profile of the central spike of 
the dark matter density and on the annihilation cross 
section of dark matter particles taking into account the 
Sommerfeld enhancement effect. 

In the presence of a small correction 5U to the New¬ 
tonian potential of the black hole, the precession angle 
of the orbit of a probe particle (SO-2 star) in one turn is 


Pir) = Ph , (2) 

where ph, rh, and /3 are the parameters. The correspond¬ 
ing total mass of the dark matter within the sphere with 
the radius r is 
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where Amin is the minimum radius to which the density 
prohle given by Eq. m expands. The subsequent calcula¬ 
tion of the precession angle of the orbit will be performed 
under the assumption that Amin < 'Pp and /3 < 3, i.e., 
that most of the mass of the dark matter within the orbit 
is located near the apocenter r = r^. We now determine 
the mass fraction of the dark matter within the orbit of 
the SO star p = [M^Mira) - MY)M{rp)]/MBK, which is 
signihcant for the subsequent analysis. 

The correction to the potential in the case of the power- 
law prohle given by Eq. m is 


(Ar^-P + ^ + C, 

\47rGp;ir^TO Inr-k^ + G, /3 = 2 


where A = 47rGp;ir^TO/[(3 — /3)(2 — /3)]. The constant G 
does not contribute to the precession angle S(j) (because 
the corresponding contribution to integral o is propor¬ 
tional to L) and the term cx 1/r is responsible only for a 
small addition to the central mass and also does not con¬ 
tribute to the precession angle. The constants Gi and 
G2 can be represented in the form Gi_2 = GmMj^yiira), 
where M-ayiiPa) is the total mass of the dark matter be¬ 
tween the event horizon of the black hole and the radius 
of the apocenter of the star under consideration. 
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The calculation of the precession angle of the orbit of 
the star in the time of one turn around the black hole 
S(j) by Eqs. m and 0 gives an expression with two con¬ 
tiguous hypergeometric functions; with the use of the 
Gauss relations for contiguous functions, this expression 
is reduced to the following expression with one hyperge¬ 
ometric function 2 ^ 1 ( 0 , b; c; z): 


6(j) 
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To test this result, we also calculated the precession an¬ 
gle 5(j) within standard perturbation theory with the use 
of the method of osculating elements [l^; the resulting 
expression coincides with Eq. 0 . The precession angle 
S(j) given by Eq. ([5]) is negative at all allowed parameters. 

The magnitude of the nonrelativistic precession angle 
given by Eq. @ is in qualitative agreement with the nu¬ 
merical calculations of precession in [l^ - [l^ . Expression 
([5|l for the precession angle at a small eccentricity of the 
orbit, e <C 1 , coincides with an accuracy of with the 
corresponding value calculated analytically by another 
method in [20|. However, at a large eccentricity, e ~ 1, 
the precession angle calculated in [ 2 ^ changes sign to 
positive and diverges in the limit e —>■ 1. The formalism 
used in is possibly applicable only at e <C 1 , be¬ 
cause the Newtonian precession angle 6(f) should always 
be negative. 

The function 5(f> given by Eq. JS]) is continuous at /3 = 2 
(see Fig. [T|). We use Eq. (I5|) to perform calculations for 
various density profiles of the dark matter. We calculate 
the function S(f>{f3,^) in Eq. ([S]) and find the level line 
6(f>{/3,^) = S(f>ohs with the value Sipohs ~ 0.01 maximum 
allowable by the observation data (see Fig.[2). The values 
P and ^ on this line indicate the parameters at which the 
observation results can be explained. 

The excess of gamma radiation from a region with a di¬ 
mension up to 10 ° from the center of the Galaxy observed 
by the Fermi-LAT telescope was fitted (see [HI, [ 2 ^) with 
the generalized NavarroFrenkWhite profile 
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where d = 20 kpc and ph(8.5 kpc) = 0.3 GeV cm“^; the 
best fit of the Fermi-LAT data is obtained at 7 ~ 1.26. 
Near the center of the Galaxy, at r d, this profile is 
close to the power-law profile 

PH(r)=po(^^^ ■ (7) 


The central region of the dark matter density peak ac¬ 
cording to Eq. o is called the cusp. If this profile is 
directly extrapolated to the center of the Galaxy, the 
mass of the dark matter within the orbit of the S2 star 
is 


Mum = M{ra) - MfR,) = 

Anpod^rl-^ - Rl-y 


2.8M0, 


( 8 ) 


where Rc = maxjr^, rann}, ?'ann being the possible inner 
edge of the distribution of the dark matter associated 
with its annihilation (see below). Quantity (jS]) is much 
smaller than the value accessible for the constraints by 
the dynamics of stars. However, in the presence of the 
central black hole, the indicated extrapolation is invalid, 
because the density profile should be significantly modi¬ 
fied by the gravity of the black hole. 

The formation of an additional density peak (spike) 
with the density profile oc r~^ around the central black 
hole was discussed in a number of works [H, [ 2 ^ . If the 
spike was formed adiabatically, i.e., gradually with an 
increase in the mass of the black hole, the density in 
the spike could be much higher than the density in the 
cusp. Fields, Shapiro and Shelton showed that /3 = 
2.36 for the adiabatically formed spike and this spike at 
{av) = const would be a very bright source at the center 
of the Galaxy (see also O). The annihilation signal from 
such a spike calculated in [ 2 ^ is a factor of about ^ 35 
stronger than the signal from the extended region with 
excess of gamma radiation. Since such bright sources at 
the center of the Galaxy are absent, the existence of the 
adiabatic spike contradicts observations. It is stated in 
[Hl| that the spike could be formed nonadiabatically or be 
destroyed. In this case, fd < 2.36 and the contradiction 
could be removed. 

Following [ 25 , we write the density of the dark mat¬ 
ter in the spike in the form of Eq. 0, where rh = 
GMualvf: ^ 1.7 pc is the radius of the action of the 
black hole, Vc = 105 ± 20 km s“^ is the observed stan¬ 
dard deviation of velocities at the distance ^ 1 pc from 
the center of the Galaxy, and the density ph is determined 
by matching Eq. Q with the density given by Eq. ([S|) at 
the radius Vh. 

The minimum radius Tann is determined by the anni¬ 
hilation of particles in the time of existence of the spike 
(see Hi). This quantity depends on the parameters of 
particles and the density distribution. According to the 
calculations in [ 2 ^ , the best fit of the gamma spectrum is 
obtained at m = 35 GeV and {av) = 1.7 x 10“^® cm^s 
in this case, rann ^ = 3 X 10“^ pc. It was shown 

in [ 2 ^ that contradiction with the bright point source is 
absent if /3 = 7 s = 1.8. In this case, the mass of the dark 
matter within the orbit of the SO-2 star is Mum — 45Mq. 
Noticeable dynamic effects (precession of the orbit of the 
star, etc.) should absent at such a small mass. 

It was assumed above that {av) = const. However, 
in a number of models of the dark matter {av) can de¬ 
pend on V. The velocities v of particles increase when 
approaching the black hole. The dependence {av) of on 
V can significantly affect annihilation. Owing to high Ke¬ 
pler velocities near the black hole, the annihilation signal 
from the center can be reduced, conserving the extended 
signal from the region ^ 10°. In particular, {av) depends 
on V in models involving the Sommerfeld enhancement 
effect (HI-[T7j. Sommerfeld enhancement is possible if a 
dark matter particle is a member of a multiplet of states 
with close masses, between which coannihilation occurs. 


3-7 
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e. g., in the model of neutralino with the dominance of 
Higgsino. The gain TZ owing to the Sommerfeld effect is 
determined from the relation {av) = TZ{av)o^ where 

7e= (9) 

Here, /r = const and b = vjc. We consider a quite general 
case where the cross section in the corresponding region 
of the parameters can be approximated by the power-law 
dependence 


(av) = {av)o , ( 10 ) 

where {crv)Q = const and Vq = const. Power-law de¬ 
pendence m was considered in in the calculation 
of the annihilation of the dark matter in self-gravitating 
bunches. The model with (av) = const and the model 
with Sommerfeld enhancement at irfi/b <IC 1 correspond 
to particular cases rj = 0 and rj = 1, respectively. 

The radius r^nn at which the maximum density p of the 
dark matter limited by the annihilation effect is reached 
is determined from the condition 

n{av)tg ~ 1, (11) 

where n = p/m and tg ~ 10 ^*^ yr is the age of the density 
peak around the black hole. For cross section m , we 
obtain 


J'ann = , 


( 12 ) 


where A = ph{crv)otg/ni. 

The velocities of particles near the black hole are 
v{r) Ri {GMBR/ry^^ at r < r/j. Let 3 — 2/3 + 7 y /2 < 0. 
This condition is satisfied for the parameters considered 
below. The corresponding rate of annihilation of the dark 
matter in the range of radii from ri to r 2 under the con¬ 
dition ri <C r 2 can be written in the form 


r2 

N = An J r'^drp^{r)m~^{<7annv) = 

m^{GMBH)'^/^{2|3-3-r]/2)' ^ ’ 

The authors of [2^, where (av) = const was accepted, 
found the parameters of the power-law profiles of the cusp 
and spike that ensure the absence of contradiction with 
the bright point source at the center. We assume that 
the integral annihilation signals from the cusp and peak 
in model m are the same as in (^ . This assumption 
allows a simple calculation without a detailed fit of the 
observed excess of gamma radiation. The main signal in 
the cusp is generated at r ~ r^, where v ^ Vc- Hence, 
fixing the parameter vn = Vc, we obtain {av)o = 1.7 x 
10“^® cm^s as in [^. The density profile in the peak 
in our case differs from that used in because of the 
dependence v(r). We determine the density profile in 



FIG. 3. Mass fraction of the dark matter ^ required by the 
observed excess of gamma radiation in the case of the anni¬ 
hilation of the dark matter particles with the cross section 
(av) oc v~'^ versus rj.). 


the peak at (av) oc v~^ taking into account the above 
assumptions. Equating the annihilation rate given by 
Eq. (IT^ in the spike at /3 = 7 s = 1.8 and {av) = (crt')o to 
corresponding rate (HI at arbitrary values [3 and {av) = 
{av)oVQ /v ’^we obtain the nonlinear equation 

^ = 0, (14) 


where 


P 


- + H + - = ^ ^ 

2 4 2 ’"^ 2js-3\f,nnJ 


(15) 


and A is determined from Eq. (US- We solve numerically 
Eq. (fT4|) with respect to x and, then, calculate the mass 
of the dark matter by Eq. ([5]). The results of the calcu¬ 
lations are shown in Eig. [3] The equality Tann = fp is 
reached at 77 = 0.6 and the inequality rann < I'a is always 
valid under the accepted conditions. Thus, in Eq. ®, 
Rc = Tp ai rj < 0.6 and Rc = rann at 7 > 0.6. The 
adiabatic density profile /3 = 2.36 is reached at 7 = 3.13. 

In particular, the mass of the dark matter Mdm within 
the orbit of the SO-2 star in particular cases 77 = 0, 1, 3.13, 
and 3.5 is 45Mo, 144Mo, 1.8 x lO^M©, and 2.8 x W^Mq, 
respectively. These values correspond to the values ^ = 
1.2 x 10-5, 3.6 X 10-5, 4.4 x 10-^, and 6.9 x 10-“^ and 
P = 1.8, 1.9, 2.36, and 2.4, respectively. These values are 
upper bounds on the possible ^ and P values. In the first 
two cases, the distributed mass of the dark matter is still 
too small to affect dynamic effects (see Eig. [3|). At the 
same time, the real prospect of the measurement of the 
additional mass of the dark matter from the precession 
of SO stars appears already at 77 > 3. 

The currently existing observation accuracy is still in¬ 
sufficient for the measurement of the precession angle of 
fast SO stars and the distributed invisible mass. How¬ 
ever, there is a high probability of reaching in the near 
future the accuracy required either for the measurement 
of the precession angle or for the determination of a 
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strong bound, which in turn will make it possible to 
impose stringent dynamic constraints on the additional 
dark mass. If the invisible mass is attributed to annihi¬ 
lating particles, the observation of the annihilation signal 
from the center of the Galaxy provides additional possi¬ 
bilities for the calculation of the distributed mass or for 
the determination of bounds on it. According to Fig. |31 
at (av) = const and even with Sommerfeld enhancement 
{av) oc 1 /u, the dynamics of stars still cannot give con¬ 
straints on annihilation, because the mass of the dark 
matter within the orbit of the SO-2 star in these cases 
is very small. At the existing accuracy, the dynamics of 


stars and annihilation are independent. However, if the 
annihilation cross section depends on the velocity with a 
large exponent 77 > 3 in Eq. m, the mass of the dark 
matter can be significant. In this case, joint constraints 
could be obtained in the near future from the dynamics 
of stars and from the data on gamma radiation from the 
center of the Galaxy. 
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